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We study a mechanism to produce the circular polarization of primordial gravitational waves. The 
circular polarization is generated during the super-inflation driven by the Gauss-Bonnet term in the 
string-inspired cosmology. The instability in the tensor mode caused by the Gauss-Bonnet term 
and the parity violation due to the gravitational Chern- Simons term are the essential ingredients of 
the mechanism. We also discuss detectability of the produced circular polarization of gravitational 
waves. It turns out that the simple model of single-field inflation contradicts CMB observations. 
' To circumvent this difficulty, we propose a two-field inflation model. In this two-field model, the 

^ , circular polarization of gravitational waves is created in the frequency range designed by the Big- 

. Bang Observer (BBO) or the deci-hertz gravitational- wave observatory (DECIGO). 

^ PACS numbers: 98.80.Cq 



I. INTRODUCTION 
. . . 

Oh' It is widely believed that the most promising candidate for an unified theory including quantum gravity is superstring 
Q \ theory. It is therefore interesting to prove or disprove superstring theory from an observational point of view. To this 
5^ ■ end, primordial gravitational waves have been considered as the most efficient probe, since the gravitational waves 

■ can carry the information of the universe at the Planck time. 

, In this paper, we focus on a robust prediction of superstring theory, namely, the parity violation due to the 
gravitational Chern-Simons term. In fact, the Chern-Simons term appears in the Green-Schwarz mechanism which is 
necessary to cancel the anomaly in the theory [l| . It also arises as a string correction . Intriguingly, the existence of 
J> ■ the Chern-Simons term can be probed by the primordial gravitational waves. This is because the Chern-Simons term 
' can generate the circular polarization of gravitational waves through the parity violation. While it is difficult to find 
OO other effects to produce the circular polarization of primordial gravitational waves. Hence, if we detect the circular 

■ polarization in the primordial gravitational waves, it would be a strong indication of existence of the Chern-Simons 
term in the early universe. Thus, it can be interpreted as an evidence of superstring theory. In addition to this 
fundamental motivation, there is a phcnomcnological interest in Chern-Simons gravity as an alternative theory of 

, gravity 0]. Even from this phcnomcnological point of view, the circular polarization of gravitational waves deserve 
detailed investigation. 

Concerning the issue if the Chern-Simons term can produce the primordial gravitational waves with circular polar- 
ization, there has been already several works It turned out that, however, there exists no observable amount 
of circular polarization of gravitational waves [7|, |8| , if the background spacetime undergoes the standard slow-roll 
inflation. The purpose of this paper is to point out that this disappointing result is reversed if we properly take into 
' account another stringy effect. The point is that the Chern-Simons term is not the only term that could be induced 
by the stringy corrections. There is another term, the so-called Gauss-Bonnet term. If we consider the Chern-Simons 
term, it would be mandatory to incorporate the Gauss-Bonnet term. Hence, we consider both terms, the Gauss-Bonnet 
and Chern-Simons terms, and study a possibility to have the circular polarization of primordial gravitational waves. 
The reason we expect a different result from the one obtained from the standard slow-roll inflation is the existence 
of the super-inflatioriary phase induced by the Gauss-Bonnet term 0, [l^, [HI, [H, [l3| ■ And more importantly, as was 
firstly discovered in [l4l . [isl . [l6j , there exists an instability in gravitational wave modes during the super- inflationary 
stage. It is this instability that generates the circular polarization of primordial gravitational waves. In fact, we show 
the primordial gravitational waves are fully polarized due to the Gauss-Bonnet term. We also discuss the detectability 
of the polarization of gravitational waves. The detectability of circular polarization of gravitational waves depends 
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on when the super-inflation occurs. Since the curvature perturbation has a blue spectrum in the super-inflationary 
regime, we need to assume the super-inflation occurs in the late stage of the slow-roll inflation. As a realization, we 
propose a two-field inflation model where the circular polarization of gravitational waves could be produced in the 
BBO 17] or DECIGO JS\ frequency range. 

The organization of this paper is as follows. In section II, we present basic equations. First, we study the evolution 
of the background spacetime in the Gauss-Bonnet-Chern-Simons gravity. For the background evolution, the Gauss- 
Bonnet term is crucial, while the Chern-Simons term plays no role. We show there are two types of inflation, namely, 
the standard slow-roll inflation and the super-inflation. Next, we derive equations for gravitational waves. Here, 
the Chern-Simons term plays an important role. In section III, we explain the essential mechanism to produce the 
circular polarization of gravitational waves. We quantize the gravitational waves in the super-inflationary background 
and calculate the degree of the circular polarization of gravitational waves. There, we will show that the significant 
circular polarization is created due to the instability caused by the Gauss-Bonnet term. We note that the model 
presented in this section is just for illustration. More realistic modes are considered in the next section. In section 
IV, we discuss the detectability of the circular polarization generated during the super-inflationary regime. We point 
out the defect of the single-field inflation model. Then, as a model consistent with current observations, we propose 
a two-field infiation model where Gauss-Bonnet and Chern-Simons term couple with the inflaton only in the second 
stage of the inflation. The final section is devoted to the conclusion. 



II. BASIC EQUATIONS IN GAUSS-BONNET-CHERN-SIMONS GRAVITY 



Before discussing primordial gravitational waves from the quantum fluctuation during the inflation, we need to 
reconsider the chaotic inflationary scenario driven by an inflaton (jj with a potential V{(j}). Although the Chern- 
Simons term does not alter the background geometry, the Gauss-Bonnet term change the inflationary scenario at the 
early stage of the inflation. The change would be important for the_primordial gravitational waves. 

We start with the action motivated from string theory given by 



S 



d'*a;^ 



-9 



d'^x^/^uj{(t))RR , (1) 



where the first term of action is the Einstein-Hilbert term and g is the determinant of the metric g^y. We have set 
the unit SttG =1. In the above action ([T]), we have taken into account the Gauss-Bonnet term 

(2) 



-^GB — R°'^'^^ RaP^S — ^R^^RaP + R^ 



and the Chern-Simons term 



RR — ^^"^''^ RappcrRjS'^ 



(3) 



where e"^^^ is the Levi-Civita tensor density. We have also allowed the coupling of the inflaton field both to the 
Gauss-Bonnet ^(0) and Chern-Simons terms ^(0). Otherwise, these topological terms vanish identically. It should 
be noted that, as is well known, the Chern-Simons term does not contribute to the dynamics of the isotropic and 
homogeneous universe. 

From now on, for simplicity, we consider a typical potential, V — l/2m^(/)^. We also specify the coupling function 
^((/)) — uj{(f>) — 160^ for concreteness, although we keep them in the formula as far as we can. For other coupling 
functions, the analysis is similar and easy to perform. 



A. Background spacetime 

Let us consider the background spacetime with spatial isotropy and homogeneity. Then, the metric is given by 

ds^ = g^^dxt'dx" = a^iri) [~dr]^ + Sijdx'dx^ . (4) 

Here, we have also assumed the flat space and used the conformal time rj. Taking the variations of the action ([1]), we 
have equations 

m' = l^'^ + lm'a^^' + ^n'e (5) 
{2n' + sn') (^1 - ^wc') +n'(i + - ^r) - m^a'^' = O (6) 

(j)" + 2H0' + 2^^^^'C,0 + m^a^<P = 0, (7) 
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where we have defined Ti. = a' /a. Here, the prime denotes the derivative with respect to the conformal time ry. 

In the conventional case, namely, without the Gauss-Bonnet term, the Hubble friction in the equation of motion 
for the scalar field makes the rolling of the scalar field slow. Then, the slow-roll inflation is soon realized for the initial 
condition with a large Hubble. However, in the presence of the Gauss-Bonnet term, the force due to the Gauss-Bonnet 
term becomes dominant for a large value of <j). If the scalar field start with a negative value, the force term accelerate 
the scalar field and makes the kinetic term dominant. Hence, the situation, Ti^ <C cj)''^ , rn^ clP ^ cf)'"^, is realized. 
Thus, Eqs.dS]) and ^ become 

a^cl)'^ + 3n^^' = (8) 

{2H' + m^)^' + H iC - n^') = . (9) 

The scalar field (f> rolls down from the negative side towards zero according to the above equations ([5]) and Now, 
it is easy to solve Eqs.® and ^ as 

cj, = -^T5/T6{-vf/\ a{r^) = {-ri)-''\ 7i = , < . (10) 

(-6r/) 

We note that the super-inflation if' > is realized in this phase. Here, H ~ H/a is the Hubble parameter. As the 
scalar field rolls down, the Gauss-Bonnet term decreases. Eventually, the conventional Hubble friction overcome the 
Gauss-Bonnet effect and the slow-roll inflation commences. A typical evolution of the spacetime is shown in Fig|T] 
In the super-inflationary phase H' > 0, the weak energy condition is violated. Hence, the system may show the 




FIG. 1: A typical evolution of the background spacetime is numerically calculated and displayed. A short period of the 
super-inflationary phase is followed by a long period of the slow-roll inflationary phase. 



instability. Of course, as you can see in Fig. 1, this instability is a transient one. To make this point clear, we rewrite 
the equations of motion ([7]) as the autonomous system for (j) and H and plot the phase flow diagram in Fig. [2l This 
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FIG. 2: Phase diagram around ~ is drawn. The thick line in the diagram denotes the trajectory corresponding to Fig. [T] 
This tells that the asymptotic solution breaks down before reaching = 0, and thus the singularity is avoided. 



diagram shows that the super-inflationary phase will be followed by the standard phase where the Hubble parameter 
is decreasing. 
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In the next section, we will show that the transient super-inflation phase gives rise to the circular polarization of 
gravitational waves. 



B. Action for gravitational waves 

In this subsection, we will deduce the quadratic action for the gravitational waves from the action ([T|). Let us 
consider the tensor perturbation 



ds^ = g^^dx^'dx" = a^{r]) [-dr/^ + (Sij + h,j)dx''dx^ 



(11) 



where hij satisfies the transverse-traceless conditions hijj — hjj = 0. Substituting the metric into the Einstein action 
with the Gauss-Bonnet term, we obtain the quadratic action jl^, [l^, [lO] 



i / d'^x^gR ~ I d'^xV^ 



^ J d^x^amlB 



2a2 



2o? 2o? 



(12) 



In addition to these terms, the parity violating Chern-Simons term gives 



1 j d^xV^u;{cP)RR^-j^ J d'xJe^'^^ 



a \b 



(13) 



where the stem | denotes a derivative with respect to the spatial coordinates. Here, we have used the convention 



Mjk 



.gijk ^123 ^ 2. It is convenient to expand hij by plane waves 



V2 



E 



d^fc 
(27r)3 



(14) 



where we have introduced the circular polarization tensor pf, {A — R, L) defined by 



„R — 



V2 



1 



(15) 



Here, the polarization tensors pf^ and p^^ are the plus and cross modes, respectively. The polarization tensors are 
normalized as p*fpE — 26^^ . The circular polarization has the property 



(16) 



where = -1-1 and A'^ = — 1 represent the right and left-handed circular polarizations, respectively. This makes the 
action of the Chern-Simons term diagonal. 

Putting above results together and defining a new variable = ' obtain 



A=R,L 



d^k 

' (27r)3 



l + ^-^jfc'lMkP + ^lMi^'^ 



(17) 



where we have defined 



A^^ = 1 . A^ 



(18) 



Thus, the equation of motion becomes 

(Mi^)" 



1 



2zi 



ZA 



(19) 



The term z'^/ za can be interpreted as the effective potential. The coefficient of the wavenumber k may be interpreted 
as the square of the speed of sound waves. 
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In the long wavelength limit, we have the solution 

„t = Gz^ + Dz{^ f^, (20) 

where G, D arc the constants of integration. In the conventional inflationary scenario, the flrst term corresponds 
to the growing mode which is actually constant when we translated back to the metric perturbations. Hence, the 
primordial gravitational waves are frozen on super-horizon scales. However, in the presence of Gauss-Bonnet and 
Chern-Simons terms, the role of the growing and decaying modes are interchanged. This gives rise to an interesting 
effect on the polarization. 



III. A MECHANISM TO PRODUCE CIRCULAR POLARIZATION 



In this section, we present a mechanism to produce the circular polarization of primordial gravitational waves. The 
interplay between the instability induced by the Gauss-Bonnet term and the parity violation due to Chern-Simons 
term is essential for the mechanism. The mechanism is efficient and generic in the sense that both terms are ubiquitous 
in string theory. 

Now, let us quantize the gravitational waves and calculate the degree of polarization. What we should do is to 
promote to the operator and expand by mode functions u^{r]) as 

H^ = a^u^ + ai\t^ , (21) 

where u^^ is the complex conjugate of and and a]^ denote the annihilation and creation operators. The mode 
functions satisfies the same equation for Once the positive frequency mode is specified, the vacuum is 
define by 

a^\0 >= . (22) 

Then, we can calculate the vacuum fluctuations 

<oiiMi^no>=iui^r (23) 

Here, we focus only on the gravitational waves which leaves the horizon during super-inflationary phase, because the 
gravitational waves still under the horizon when super-inflation ends are the same as th ose of the slow-roll i nflation. 

During the stage prior to the slow-roll inflation, za can be approximated as za = y'-ne/^- X^kuj'/2. In spite 
of this simplification, it is still difficult to solve the equation of motion analytically. Hence, we solve the equation at 
super-horizon scale and sub-horizon scale, separately. After that, we connect these solutions smoothly at the horizon 
crossing. 

First, let us consider the sub-horizon scale, H <^ k ov —krj ^ 1/6. Leaving the terms up to the order of l/(— fcry), 
we obtain the equation 

(u^y' + k^^l-X^l^y^ = 0. (24) 

This can be solved analytically by using confluent hypergeometric functions. Instead of doing that, we estimate the 
degree of the polarization approximately. The initial condition can be set by imposing the Bunch-Davies vacuum. 
More precisely, we choose the positive frequency mode 



in the asymptotic past. Namely, deep inside the Hubble horizon, both the left and right polarization modes are simply 
oscillating. At the time —kr] = 8/3, the asymmetry shows up. After that time, —kr] < 8/3, we can approximate the 
solutions as 



<=^^^^k^ ^^M^V— (-^^^r^H^P^ exp -2^— (-M , (26) 
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where we have used WKB approximation which is vahd for the period we are considering. To determine the constants 
of integration Ai and A2, we have to match solutions smoothly at —k-q^ 8/3 as 



J)<8/3 



<(^)|,>8/3' «) (^)|,<8/3 = «) (^)|,>8/3' 



(27) 



Apparently, both Ai and A2 are not zero. Hence, for the left-handed circular polarization modes = —1, the 
solution shows damping oscillation. While, for the right-handed circular polarization mode = 1, the solution 
grows rapidly. It is this instability due to the Gauss-Bonnet term that produces the difference between the left and 
right-handed circular polarization modes. The ratio can be quantified by 



|ur|2 exp(-8/3) 



exp(-32/3) 



2980 



(28) 



Thus, we can expect the fully polarized gravitational waves. 

To make the analysis more accurate, we need the solutions to continue to the super-horizon scale, Ti ^ k or 
— krj <^ 1/6. On the super-horizon scales, the solution can be obtained as 



0. 



This can be solved easily as 



Bi{-kf])^/^ + B2{-k?]) 



2/3 



(29) 



(30) 



where Bi and B2 are constants of integration. To determine the constants Bi and B2, we have to match solutions 
smoothly at — r/cross = l/(6fc) as 



"k('7cross)|.^^j, — U^{Vcross)\-^^i. , ( ""k )' ('ycross) | .^^j, — ( Wj^ )' (?7cross ) | -^^j. • (31) 

What we want to calculate is the degree of the circular polarization defined as the difference between the power of 
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FIG. 3: The degree of the polarization n(fe) as a function of wave numbers k is shown. As expected, the gravitational waves 
are almost 100% circularly polarized. 



right and left-handed circularly polarized gravitational waves at the end of the super inflation phase: 



n(fc) 



Kfand)p-k^(r;cnd)P 

K(r?c„d)P + |u{^(r/cnd)p 



(32) 



where jycnd represents the time when the super-inflation ends. We have performed matching analysis and displayed 
a numerical result in Fig[31 Because of the behavior just before the horizon crossing, namely right-handed circular 
polarization mode is growing and left-handed circular polarization mode is decaying, consequently the net polarization 
is produced. The resultant polarization is sufficiently large and hence detectable by BBO or DECIGO. 

To calculate the spectrum of the polarization of gravitational waves which is directly observed today, we have 
to connect solutions from super inflation phase to the slow- roll inflation, radiation dominant, and matter dominant 
phase. This part of calculations is standard and does not alter the polarization. 
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Finally, we should mention a subtle point in the calculation. From the definition of za (jlSp . it is obvious that there 
is a point rising where = 0. That occurs close to the time when the left-handed circular polarization mode cross 
the horizon. There, the solution behaves as 



"k = C'lyh-'fengl +C2y h-?7sing|ln(|?7-?7sing|) , (33) 

where Ci and C2 are constants of integration. We can see that 0, in the limit rj T?sing, so it causes no problem 

to solve the Eg. pop . However, the relevant quantity is the physical amplitude ■0^, rather than namely, 

^k = — -Cl+C2H\V~Vsins\)- (34) 

Apparently, the physical amplitude diverges at 77 ^ '7sing- In the previous work [71, this issue was serious, because 
this divergence means the breakdown of the linear analysis. It suggests that we should take into account the non- 
linear effect which is not available. In our case, the relevant mode is the right polarization mode which has no 
divergence problem. As the breakdown of the linear analysis occurs just at the horizon crossing, the non-linear effects 
works almost instantaneously. Hence, it is reasonable to assume that the non-linear effect for the left-handed circular 
polarization mode hardly affects the evolution of the right-handed circular polarization mode. Therefore, we believe 
the consequence about the degree of polarization does not change. 



IV. DETECTABILITY OF CIRCULAR POLARIZATION OF GRAVITATIONAL WAVES 



First of all, we should emphasize that the analysis done in the previous section is explanatory. The frequency range 
where the circular polarization is created will change depending on when the super-inflation occurs. In this section, 
we will discuss the indirect detection through CMB and the direct detections by interferometer. 

As to the in-direct detection of circular polarization through the cosmic microwave background radiation (CMB), 
the required degree of circular polarization has been obtained as |H| > 0.35(r/0.05)~°-^ [21|, where r is the tensor- 
to-scalar ratio. The relevant frequency in this case is around / ~ 10""'^'^ Hz. As to the direct detection of circular 
polarization, the required degree of circular polarization has been estimated as H ~ 0.08(ilGw/10^^^)^^(SNR/5) 
around the frequency / 1 Hz [l^, [^^l , where flow is the density parameter of the stochastic gravitational waves 
and SNR is the signal to the noise ratio. Here, 10 years observational time is assumed. Of course, for the isotropic 
component of circular polarization of stochastic gravitational waves, one needs to break the symmetry of the detector 
configuration. 
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FIG. 4: Tlie numerical result of the spectral index n for the power spectrum of primordial gravitational waves as a function of 
wavenumbers k is shown. In the relevant frequency range, we have obtained the spectral index n — 2.66 numerically. 

Thus, in order to know if we can detect polarized gravitational waves, we need the degree of the polarization, the 
intensity of the gravitational waves J^gw or the tensor-to-scalar ratio r, and also the spectrum around the relevant 
frequency ranges. As we have shown already, we can get the polarization H ~ 1 easily. Moreover, it is not difficult to 
get the required fiow or r provided standard model parameters. We can calculate the spectrum numerically assuming 
the Bunch-Davis vacuum. We have presented a numerical result in FigU) This result can be understood analytically. 
In fact, the initial mode function behaves oc 1/Vk inside the horizon which should match to the growing mode 
Bi{—kr]y^^ on the super horizon scales. Matching condition at —kr] ~ 1 gives the k dependence of the coefficient Bi 
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as Bi cx 1/Vk. Thus, we have Uk oc k^/^-'^/^ ^ Consequently, the power spectrum for the gravitational waves 

becomes /c'^Im^P oc fc^/^. Since the spectrum has the blue index n ~ 2.66, the polarization would be detectable even 
in the low intensity case. Now, the point is the frequency range, namely, the time when the polarization is created. 
In the following discussion, we will mainly focus on this issue. 

In the case of the single inflaton model we have discussed, the frequency where the circular polarization is significant 
lies around the CMB scale / ^ 10^^^ Hz. Hence, there is a possibility for the circular polarization to be indirectly 
detected through the CMB observation, more precisely, through the Temperature and B mode or E mode and B 
mode correlations [1, [2l|, In a recent paper [2l|, it is argued that a high degree of polarization is necessary 
for the polarization to be detected in the future observations. The high degree of polarization can be attainable by 
the mechanism we have found in this paper. However, in this simple model, we have to worry about the curvature 
perturbations which is also expected to have the blue spectrum during the super-inflation. One possible way out of 
this problem is to assume the curvature perturbation created by the inflaton is negligible and to resort to the curvaton 
scenario [1^ [H, as is common to the pre-big-bang type models [H, [2^ [sol l3ll|. However, it is not clear if the 
concrete realization of this kind of model is possible. Moreover, one may feel that introducing the curvaton mechanism 
into the scenario is not appealing even if it exists. Fortunately, there exists another interesting possibility. We will 
devote the remaining part of this section to the model which is consistent with observational constraints. 

To circumvent the problem of the curvature perturbations, we propose a two-field inflation model defined by the 
action 



S = ^ I dli 




1 

16 



A^xy/^uj((l))RR , (35) 



where the potential V{<j)^x) is assumed to take a hybrid type. The main difference between our model and the 
conventional inflation model is simply the Gauss-Bonnet and Chern-Simons term. In the first stage, the field x drives 
the slow-roll inflation. During this stage, the field takes a constant large field value. Hence, both the Gauss- 
Bonnet and the Chern-Simons terms are decoupled from the system. During this conventional inflationary period, 
the temperature fluctuations of CMB are created. And the spectrum of the fluctuations is almost flat. Then, at 
some point, the motion in the direction of </> is triggered and the Gauss-Bonnet and Chern-Simons terms take part 
in the dynamics. In this second stage, the (j) field rolls down the hill rapidly and the Gauss-Bonnet term induces a 
super-inflationary phase. During this stage, the circular polarization of the primordial gravitational waves is produced 
through the parity violation of the Chern-Simons term. Subsequently, the second standard slow-roll inflation driven 
by the same field takes over. 

Here, we present an example. The potential function is 




FIG. 5: Potential function for the two field model is depicted. Initially, the scalar field (ji is sticked to the point 4> = —3000 
and the other scalar field x slowly rolls down to give rise to the infiation. At some point, <j) starts to roll down toward <j) = 0. 
Then, the coupling to the Gauss-Bonnet term induces the super-infiation where the circular polarization of gravitational waves 
is created. Subsequently, the standard slow-roll infiation foUows. 
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y (</>, x) = {\x^ + 1) { (2/(/.,„)2 [(</> + (/.„V2)2 - (</.„72)2] } 




where the parameters are set as 

(^m = 3000, = 10"^ = 0.235, = 10^ = lO"*' . (37) 

This potential is shown in Fig. [S] Couphng ^ and w are same as those of the single- field model, namely ^ = w = 16(/)'*. 
We have calculated the evolution of Hubble parameter numerically (see Fig. |6l). 
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FIG. 6: The time evolution of Hubble parameter is numerically calculated. The curvature perturbations on the CMB scales 
are produced during the slow-roll phase. Hence, this model is consistent with CMB observations. Gauss-Bonnet term induces 
the super-inflation around TV ~ 50 where the circular polarization would be created. It could be observed directly through the 
interferometer detector. 
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FIG. 7; A schematic picture of the expected spectrum is depicted. The frequency range where we can observe the circular 
polarization depends on the model parameters. In that range, the spectrum is blue. In other parts, the spectrum is almost 
flat. The expected sensitivities of BBO and DECIGO are also plotted. 

For this model, the duration of the first and second slow-roll infiation is tunable. Hence, it is possible to have the 
gravitational waves with the significant circular polarization in the BBO or DECIGO frequency range (see Fig[7]). 
Since the spectrum is blue, even the ground based detectors may find the circular polarization as is pointed out 
recently [32l |. 

To make a more precise prediction of the polarization of gravitational waves, we need more complete analysis. In 
particular, we should calculate the spectrum of the curvature perturbations of this model. We leave these details for 
future publications |33j]. 
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V. CONCLUSION 



We have studied a mechanism to produce the circular polarization of gravitational waves in the string-inspired 
cosmology. It turned out that the circularly polarized gravitational waves are ubiquitous in string cosmology. The 
point we have observed is that there are two key terms in string theory, namely, the Chern-Simons term and the 
Gauss-Bonnet term. The Chern-Simons term violates the parity invariance, therefore it makes a room for the circularly 
polarized gravitational waves to be produced. However, in the previous works [3, 0)11] j it had been shown that there 
is no significant circular polarization of gravitational waves within the conventional inflationary scenario. In this 
paper, we have shown that the Gauss-Bonnet term reversed the previous conclusion. The Gauss-Bonnet term has 
changed the background evolution in such a way that the super-inflationary epoch appears during the conventional 
inflationary stage. During the super-inflation, there exists an instability in the tensor modes. It is the instability that 
produces a significant circular polarization 11 ^ 1. We have also discussed the detectability of the polarization. We 
have shown that the detectability depends on the specific scenario of inflation. In the single-field inflation model, the 
circular polarization would be observable at the CMB scale. The consistency with the observation of the temperature 
fluctuations requires the curvaton scenario in this case. To circumvent the situation, we have also proposed a two-field 
inflation model. In the two-field model, the frequency range where the circular polarization is produced depends on 
the parameter. It would be possible to make a model to produce gravitational waves with circular polarization in the 
BBO or DECIGO frequency range. In this case, curvature fluctuations on large scales are produced by the standard 
slow-roll inflation. Hence, we do not need the curvaton scenario any more. 

There are several directions to be explored. We should investigate the two-field inflation model in detail [s^]. In 
particular, it is interesting to seek a concrete realization of the model in the string landscape. It is also intriguing to 
search other mechanism to produce a circular polarization of gravitational waves. One simple possibility is to consider 
the non- vacuum initial state in Lorentz violating inflationary scenario [s^]. The coupling to the Lorentz violating 
sector could violate the parity invariance. Hence, it is possible to take the asymmetric vacuum for right and left-handed 
circular polarizations. In fact, the Chern-Simons term already allows this possibility. There may be other interesting 
mechanisms with non-inflationary origin, such as the turbulence or the helical magnetic field s |24 | . |36|. The Gauss- 
Bonnet term at present may become an interesting subject in the light of the dark energy [stI. IssI. [39l. |40|. l4ll. |4^. l43l |43| . 
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